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Abstract. X-ray observations of SN 1987A in the Magel- 
lanic Clouds have been performed throughout the ROSAT 
mission, using the PSPC and the HRI detectors. We 
present the X-ray light curve based on all observations 
in the years 1991-1995. For the first time a significant in- 
crease of the X-ray flux from SN 1987A can be seen, cor- 
responding to X-ray luminosities (0.5-2 keV) of 0.8 — 2.2 x 
10 34 erg/s about 4-8 years after the explosion. SN 1987A 
is surrounded by a ringlike nebula, which is thought to 
form the interface between the blue-supergiant wind and 
the denser red-supergiant wind of the progenitor. The X- 
ray data can constrain the density of the matter inside 
the ring to about 30 amu cm~ 3 and the date at which the 
blast wave will reach the ring to about AD 2003, when a 
dramatic brightening is expected to occur. Nevertheless, 
other interpretations of the X-ray emission are possible. 
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stars - supernova remnants 



1. Introduction 

The explosion of SN 1987A (Shelton 1987) was one of 
the exciting historical astronomical events. The neutrino 
flux detected at the time of the explosion (Koshiba et al. 
1987), and the observations of X- and soft 7-rays (Sun- 
yaev et al. 1987, Dotani et al. 1987) constitute some of the 
"first observations" of such phenomena from supernovae 
confirming the general correctness of our theoretical un- 
derstanding. In addition SN 1987A is classified as type II, 
therefore a neutron star is expected at the center of this 
remnant. 

Hard X-rays have been detected from SN 1987A with 
the Kvant module on the Russian MIR station and with 
the Japanese Ginga satellite (Sunyaev et al. 1987; Dotani 
et al. 1987). The high-energy X-rays originate from the 
decay of 56 Co, where the 7-ray lines undergo multiple 
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Compton scatterings and are degraded into X-rays which, 
however, are absorbed through the tenuous cocoon of the 
stellar ejecta which is thick enough to block out soft X- 
rays from the compact object for quite a long time. 

The "soft" X-ray component originally observed with 
Ginga below ~ 10 keV could not be confirmed by 
MIR-TTM and MPE rocket observations (Aschenbach 
et al. 1987), nor by the ROSAT first-light observations 
(Triimper et al., 1991). Since 1991, however, the super- 
nova is detected as a faint soft X-ray source (Beuermann 
et al., 1994; Gorenstein et al., 1994). In this Letter we 
present data from continuous monitoring of SN 1987A 
with ROSAT. 



Table 1. Observation Summary 



Date 


Day 1 


Instr. 


Time 


Count Rate 2 








[s] 


[10" 3 c/s] 


16.06.90-28.06.90 


1215 


PSPC 


6400 


< 2.7 a 


11.02.91-13.02.91 


1448 


HRI 


23465 


0.8 ± 1.0 


21.04.91-07.10.91 


1645 


PSPC 


24908 


1.8 ±0.5 


03.02.92-14.05.92 


1872 


PSPC 


27288 


2.2 ± 0.5 


04.12.92-05.07.93 


2258 


PSPC 


23404 


3.2 ±0.6 


28.09.93-30.09.93 


2408 


PSPC 


9426 


3.8 ±0.8 


23.12.93-04.10.94 


2715 


HRI 


21182 


3.9 ± 1.3 


01.01.95-11.10.95 


3013 


HRI 


62420 


4.8 ±0.7 



mean day after explosion 

2 PSPC channel 52-201, HRI channel 2-8 multiplied with 2.65 

3 from Triimper et al., 1991 



2. Observations and Results 

SN 1987A was the target of the ROSAT first light observa- 
tion, and was scanned during the ROSAT All-Sky Survey. 
In both cases only upper limits could be derived (Triimper 
et al., 1991). Since then SN 1987A has been monitored reg- 
ularly with both the PSPC and the HRI detectors aboard 
ROSAT. Table 1 gives an overview of the observations. 

For the purpose of this paper all ROSAT data of 
SN 1987A obtained since the beginning of 1991 have 
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Fig. 1. X-ray image of the region of SN 1987A with the PSPC 
(top) in the pulse height range 52-201 and the HRI (bottom) 
in the pulse height range 2-8. 



been analysed using the MPE interactive analysis system 
EXSAS (Zimmermann et al., 1994). For the PSPC data 
all good time intervals selected by the ROSAT standard 
analysis (SASS; Voges et al. 1992) have been analysed, 
yielding a net observing time of 85026 seconds. HRI events 
have been selected using a custom made procedure after 
having joined the SASS accepted and rejected event files. 
This resulted in a net observing time of 107066 seconds, 
a gain of ~ 5% compared to the SASS products. Figure 
1 shows the X-ray images centered on SN 1987A. The 
PSPC image has been accumulated in the energy range ~ 
0.5-2 keV (pulse height channels 52-201). The HRI image 
was restricted to the pulse height interval 2-8, accepting 
a minimal loss (~ 6%) of X-rays but at the same time sig- 



nificantly reducing the intrinsic and particle-induced HRI 
background (see David et al., 1995). The supernova shows 
up as a weak point source embedded in a rather complex 
region of diffuse X-ray emission. The average HRI position 
i?. A (2000) = 5' l 35 m 28.5", £>£:C(2000) = -69°16'11.5" is 
within 2 arcsec from the catalogued position; no evidence 
is seen of a possible second, nearby X-ray source suggested 
by Gorenstein et al. (1994). 



Table 2. Spectral Analysis 



Model 1 




dof 


r/kT 

[keV] 


Flux 2 / EM a 


power law 


0.65 


5 


3.3 ± 0.4 


5.38 ±0.81 


thermal brems 


0.96 


5 


0.44 ±0.11 


5.36 ±0.84 


blackbody 


1.74 


5 


0.17 ±0.02 


7.88 ±2.0 


RS (z=l) 


3.61 


5 


1.00 ±0.16 


(5.0 ± 1.5) • 10 56 


RS (z=0.3) 


3.00 


5 


0.99 ±0.17 


(1.4 ±0.4) • 10 57 


RS (z=3) 


3.87 


5 


1.00 ±0.17 


(1.8 ±0.6) • 10 56 


2T-RS (z=l) 


1.22 


3 


0.14 


1.8 ■ 10 57 








1.28 


6.5 ■ 10 56 



1 N H fixed to 1.1 

2 Flux in 10~ 14 



erg cm 
Emission measure in cm" 



1 (PL, TB and BB) 
(RS models) 



> 

CD 



C 
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Fig. 2. PSPC pulse height spectrum of SN 1987A. The solid 
line shows the model for a 0.44 keV thermal bremsstrahlung 
spectrum and a neutral hydrogen column density fixed to 



For the spectral analysis care has been taken of the 
complex diffuse emission in the vicinity of the supernova. 
PSPC source photons were accumulated in a circle around 
SN 1987A with radius 36" , while background photons were 
taken from an annulus with inner radius 36" and outer ra- 
dius 82" . The pulse height distribution in the range 11-250 
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channels has been binned to a signal to noise ratio of 4, 
yielding 318 net counts. Various spectral models have been 
fit to the data, assuming a fixed line-of-sight absorption 
column density of Nh — 10 21 cm" 2 (Gorenstein et al., 
1994). The PSPC response matrix is known to change as a 
function of time (Prieto et al., 1996), a fact which is taken 
care of by EXSAS using different instrument matrices for 
the beginning and the end of the PSPC lifetime. Since the 
PSPC data covers a long period of time (1991 Apr - 1993 
Sep) the fits were performed with both response matrices 
in order to test for systematic errors, however no signifi- 
cant differences were found. Table 2 gives an overview of 
the spectral fits. Continuum spectra like a simple power 
law with photon index T = 3.3 or thermal bremsstrahlung 
with kT — 0.44 keV fit the data well. Figure 2 shows 
the PSPC pulse height spectrum with the corresponding 
thermal bremsstrahlung model. The average 1 0.5-2 kcV lu- 
minosity derived for the power law or bremsstrahlung fit 
is (1.6 ± 0.2) • 10 34 erg s _1 , assuming a distance of 51 
kpc to the LMC. A single-temperature Raymond-Smith 
thermal plasma model (RS), which is dominated by few 
emission lines, is rejected at a probability greater than 
99%, even if the abundance is varied. This may indicate 
a more complicated ionization or temperature distribu- 
tion. A two-temperature Raymond-Smith model is given 
for completeness, however, the errors are too large to draw 
firm conclusions. 



In order to obtain a light curve of SN 1987A, sepa- 
rate images were produced for all time intervals in table 
1. Source counts were extracted from a ring around the 
supernova with radius 24" for the HRI and 40" for the 
PSPC, respectively, containing roughly 90% of the HRI 
(David et al., 1995) and about 95% of the PSPC photons 
(Hasinger et al., 1992), respectively. Background counts 
were extracted from an annulus around the source ring, 
having an outer radius of 82". Net source count rates 
were corrected for deadtime and the point-spread-function 
losses. In order to compare the HRI and PSPC count rates 
we computed the ratio of the count-to-flux conversion fac- 
tors for both instruments, taking into account the selected 
pulse height channel range (channel 52-201 for the PSPC 
and channel 2-8 for the HRI). Assuming the first three 
spectral models in table 2, the predicted count rate ra- 
tio factor varies between 2.46 and 2.89. We therefore as- 
sumed a PSPC/HRI conversion factor of 2.65. The final 
equivalent PSPC count rates are given in table 1. These 
count rates can be converted to 0.5-2 keV fluxes assum- 
ing a factor of 6.95 • 10 10 PSPC counts per erg cmT 2 s" 1 . 
Figure 3 shows the X-ray light curve of SN 1987A, clearly 
indicating a continuous flux increase during the interval 
1500-3000 days after the explosion. We fitted the count 
rates with a power law model as a function of time and 
derive an exponent of 1.67 ± 0.35 (lc). A linear increse 
starting on day ~ 900 can fit the data as well. 



ROSAT observations of SN1987A 
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Fig. 3. Time evolution of the X-ray emission of SN 1987A 
The data points marked with filled circles correspond to PSPC 
observations (channel 52-201), open circles are from HRI ob- 
servations (channel 3-8) multiplied by a factor of 2.65. The 
dashed line indicates a quadratic fit to the data. A linear in- 
crease starting on day ~ 900 can fit the data as well (dotted 
line) . 



3. Discussion 

Using the FOC aboard the Hubble Space Telescope a faint 
circumstellar nebula has been observed around SN 1987A 
which has the shape of a thin ring with a radius of 0.2 
pc and a thickness of only 0.02 pc (Jakobsen et al. 1991). 
The authors estimate a density of n e ~ 2 • 10 4 cto~ 3 and a 
total mass of M <~ 0.2M Q for the ring, which is thought 
to originate from a shock front in the progenitor wind at 
the interface between a slow, tenuous red supergiant wind 
(RSW) and a fast asymmetrical low-density blue super- 
giant wind (BSW). In the meantime high quality images of 
the Hubble Space Telescope show two fainter rings above 
and below the original ring nebula which probably mark 
the hour-glass shaped transition region between the BSG 
and the RSG winds, where a substantial density enhance- 
ment is expected to occur. The ring nebula is also observed 
to expand very slowly, presumably somewhat faster than 
the RSG wind speed (Plait et al., 1995). 

Several investigators have discussed the X-ray emis- 
sion of supernovae, in particular that of SN 1987, in terms 
of thermal radiation of a reverse shock wave originating 
from the interaction between the supernova ejecta and the 
previously existing circumstellar material (e.g. Chevalier 
1982). The supernova blast wave is expected to interact 
with the wind, first with the BSG and later with the RSG 
wind, causing a shock wave to heat up the circumstellar 
material to X-ray temperatures. 
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The X-ray luminosity of the shock wave and its tem- 
poral behaviour depend on the density and gradient of the 
ambient medium. In the case of a freely expanding BSG 
wind with a density p ~ r~ 2 a decline of the X-ray lu- 
minosity as a function of time according to Lx ~ t~ ' 87 
is expected (Beuermann et al., 1994; Gorenstein et al., 
1994). Beuermann et al. have already remarked that the 
measured X-ray luminosity is too high for the simple free 
wind model. Our current X-ray light curve is clearly in- 
consistent with the free wind model. 

However, there are reasons to believe that a simple 
r~ 2 density dependence is not the correct description for 
the BSG wind of SN 1987A. As the existence of the ring 
nebula and the terminating discontinuity prove, the BSG 
wind must catch up and interact with the RSG wind and 
is therefore no longer freely expanding. Several authors 
recently assumed the BSG wind matter to be homoge- 
neously distributed within the radius of the ring nebula 
(Suzuki et al., 1993; Masai and Nomoto; 1994; Luo et al., 
1994). In this case a steadily increasing X-ray luminosity is 
predicted in the BSG wind interaction phase, with a time 
dependence proportional to t 2 04 , which is consistent with 
our findings (see fig. 3). Chevalier & Dwarkadas (1995) 
propose that the supernova shock is now interacting with 
an Hu region created by the BSG star in the swept-up 
RSG wind. They predict a slow rise of the X-ray flux, 
which, however should start only around day ~ 1000. Our 
supernova X-ray light curve is consistent with this model, 
too. 

Masai & Nomoto (1994) present detailed hydrodynam- 
ical calculations and predict a dramatic increase in X-ray 
luminosity around the time when the stellar ejecta hit the 
dense ring nebula. Because the SN blast wave is slowed 
down by the shock interaction, the impact time depends 
on the density of the ambient medium according to t r i ng ~ 
p 1 / 6 . Different authors assume substantially different val- 
ues for this density, ranging from p w 2 amu cm~ z (Luo et 
al., 1994; Masai & Nomoto, 1994) to p w 100 amu cm~ 3 
in the equatorial plane (Chevalier & Dwarkadas, 1995). 
Correspondingly, the range of predicted impact dates is 
still considerable (roughly 1996 - 2008). The new ROSAT 
determination of the X-ray lightcurve of SN 1987A now 
allows to put further constraints on the density in the cir- 
cumstellar material. The total X-ray flux measured on day 
3000 (see fig. 3 and tab. 1) corresponds to a 0.5-2 keV lu- 
minosity of L( . 5-2) = 2.0-10 34 erg s _1 . Depending on the 
assumed spectral model, the energy band correction factor 
is quite uncertain. Assuming a correction factor of 3 for 
a thermal bremsstrahlung spectrum we estimate a total 
X-ray luminosity of Lx = 6 • 10 34 erg s _1 . If we interpret 
this as coming from the SN blast wave in the context of 
the Masai & Nomoto (1994) model, we can estimate the 
density of the BSG wind to pbsw ~ 30 amu cm~ 3 from 
their equation (7). Correspondingly, the impact time of 
the SN blast wave at the location of the ring nebula would 
be t r ing ~ 15.7 yr after the explosion, i.e. in the year 



AD 2003. We have to caution, however, that we find sig- 
nificantly lower temperatures than predicted by Masai & 
Nomoto. A factor of ~ 3 higher density and, correspond- 
ingly, a later impact time has been proposed by Chevalier 
and Dwarkadas (1995). In view of the rather complicated 
situation these predictions are still quite uncertain. 

There could be other contributions to the observed X- 
ray emission of SN 1987A, in particular radiation from 
the surface or the magnctosphcrc of a putative neutron 
star. Actually, the currently observed X-ray luminosity 
and temperature are comparable to the values expected 
from the thermal surface radiation of a young neutron 
star. However, the circumstcllar cocoon of matter is ex- 
pected to block out all but the highest energy X-rays of 
the neutron star for a long time. Under relatively opti- 
mistic assumptions, Kumagai et al. (1993) predict a pos- 
sible pulsar to be seen after ~ 8 yr only at energies above 
6 keV. Nevertheless, relatively little is known about the 
dumpiness of the absorber and thus the observed increase 
of X-ray flux could still be associated to radiation from 
a compact object. The long-term X-ray lightcurve of SN 
1987A can likely discriminate between the different pos- 
sibilities. Continuing X-ray observations of SN 1987A are 
therefore highly desirable, in particular since at any rate 
a major brightening is expected in the coming years. 
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